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INTRODUCTION 

Given a polynomial vector field of the plane: 



X = y)d/dx + Q{x, y)d/dy, 

where P and Q are both polynomials of degree less than d, Hilbert's 16th problem, part 
B, calls for finding a bound (which depends only on d) to the number of limit cycles 
(isolated periodic solutions) of X . Despite the many contributions to this question, the 
problem still remains open nowadays even in the case d = 2. One of the most spectacular 
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contribution to this field is due to N.N. Bautin ([1]). These notes aim at providing 
generalizations of Bautin's theorem which involve several new ideas introduced in recent 
years. The creation and developments of the theory of non linear oscillations is due to the 
school of A.A.Andronov. The famous books " Theory of oscillations " by Andronov, 
Khaikin, Witt" followed by "Qualitative theory of second order dynamical systems" 
and "Theory of bifurcations of dynamical systems in the plane" run several editions in 
various languages. As one of the most prominent followers of Andronov, N.N. Bautin 
wrote in 1939 a masterpiece contribution to the classical center-focus problem posed by 
Poincarc. In his lecture given on the occasion of Arnol'd's 60-th birthday at the 
Fields Institute in Toronto, June 1997, ([13], [14]) S. Smale listed 18 problems chosen 
with these criteria: simple statement, also preferably mathematically precise; personal 
acquaintance with the problem ; a belief that the question, its solution, partial results, or 
even attempts at its solution are likely to have great importance for Mathematics and its 
developments in the next century. Besides the Riemann Hypothesis, one was on Hilbert's 
list; the Hilbert's 16-th problem on limit cycles. The author hopes that the specialists 
will be interested to have here collected techniques which appeared in a scattered way 
(from 1995 to 2000) and have been revisited in a more systematic manner. For non- 
specialists, this survey may contribute to introduce this vivid subject which involves 
so different methods such as projections of analytic sets, division by an ideal, complex 
foliations, Riemann surfaces, algebraic invariant theory and normal forms, and which 
belongs to the tradition of the qualitative theory of differential equations. Chapters I-II 
review results which have essentially already been published jointly with Yosi Yomdin 
([8], [9]). Chapter III is perhaps the most original and develops an extension of Bautin's 
approach to any dimension in the framework of the theory of pertubations of a Nambu 
dynamics ([12]). 



2 



I- CLASSICAL BAUTIN's APPROACH 



The first paragraph is devoted to the presentation of Bautin approach ([1]). We 
give the definitions of the center set, of the Bautin ideal and of the Bautin index. We 
show that the center set is a real algebraic variety. We select particular generators of 
the Bautin ideal which arc adapted to the first return map in a sense to be precised. 
We give an estimate of the domain of convergence of the return mapping and we show 
accordingly to ([2], [3], [4]) that it is an ^o-Series. We conclude this part with Bautin's 
argument based on the classical RoUe property to produce a bound to the number of 
real zeros of the return mapping. In this survey, we will mostly be concerned with 
homogeneous nonlinearities. In this paragraph, it is not necessary to assume that the 
perturbation is homogeneous: all definitions and theorems can be easily modified with- 
out any fundamental changes. Consideration of this particular case allows simplification 
of the notations and also it captures the essential properties of the general case. 

I.l Construction of the return mapping by Bautin's method 

We consider a polynomial plane vector field X of type: 

X = xd/dy — yd/dx + ^ [aijx'^y^ d/dx + bijx'^y^ d/dy]. (1.1) 

i,j/i+j=d 

and the associated fiow, solution of the system: 

x=-y+ ^ ttijx'y^ = -y + P{x,y), (1.2a) 

i,j/i+j=d 

y = x+ hj^y = -y + Q{x, y). (1.26) 

The parameters of the vector field (a, h) can take any real values and thus (a, h) 
should be considered as a point in the vector space B?'^'^^^\ We first recall Bautin's 
approach to find the return mapping of (1.1) near the origin. Write (1.2) in polar 
coordinates (r, d): 

3 



X = rcos(^), y — rsm{9). 



(1.3) 



This leads to: 

2rr = 2{xx + yy),rr = xP + yQ = r'^+^A{9), (1.4) 

e^{xy- yx)/{x^ + y^) = l + r'^-^B{e), (1.5) 

where A{9) and B{9) are two trigonometric polynomials (in cos(^), sin(^)) linear in 
the parameters (a, h). 
This yields: 

dr/dd = r'^A{e)/[l + r'^-^B{e)], (1.6) 

and thus 

oo 

dr/de = ^{-lfr^^'^-^^+'^A{e)B{ef. (1.7) 

A;=0 

This equation (1.7) may be rewritten: 

dr/de = ^r^Rk{e), (1.8) 

k>d 

where the coeflBcients Rk{0) are trigonometric polynomials in (cos(^), sin(^)) and 
polynomials in the parameters (a, h). To simplify the notations, the dependence on the 
parameters (a, h) is not made explicit. Bautin looks for a solution of (1.8) r = r{9) so 
that r(0) = ro, given as an expansion: 



r = ro + V2{e)rl + ... + VkiOy'^ + . 

Comparison between (1.8) and (1.9) yields: 
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(1.9) 



(1.10) 



v'M = Rd{9), 



(1.11) 



k 



ae) = J2Bik[vdie),...,vk-iie)]Riie),k>d + i. 



(1.12) 



The polynomial Bij.[ad, ak-i] displays integer coefficients and is the coefficient of 
X'^-* in {X+adX^^+.-.+apXP+.-.y. The relation (1.12) allows to determine inductively 
the functions Vk{d)'- 



Two facts can be easily derived from this construction: 

i) Vk{0) is polynomial in 9 (of degree less than k) and in {sin{9) , cos{9)) . 

ii) Vk{9) is polynomial in the parameters (a, b) of the vector field. Thus in particular 
the coefficients ffc(27r) of the return mapping are polynomials in (a, b). 

Definition 1. 1.1 

The vector field X displays a center at the origin (0 e i?^) (or is said to be a center) 
if and only if r(27r) = r(0) for all values of r(0) close to 0. 

1.2 Convergence of the Taylor series, majorant series and Ao-series 

Let fx{x) = Y^o^kWx^ be an analytic series in x with polynomial coefficients in 
the parameters A — (Aq, Ad). Denote | | (norm of the polynomial afc(A)) as the 
sum of the absolute value of the coefficients. 

Definition 1.2.1 

The series fx is called an ylo-series if the following two conditions are satisfied: 
There are positive constants ^1,^2,^3, K4 such that: 




(1.13) 
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1- deg{ak) < Kik + K2, 

2- I ak \<K3Kl 

Ao-series form a subring of the ring of formal power series in x with polynomial 
coefficients in A. All the usual analytic operations, like substitution to a given analytic 
function, composition, inversion, etc... transform Ao-series into themselves. (The proof 
is rather straightforward and is not provided here). The ^o-series have been precisely 
introduced (in the subject) by M. Briskin and Y. Yomdin ([4]). The proof of the 
following easy lemma is also left to the reader: 

Lemma 1.2.2 

An ^o-series f\{x) converges in the disc D(0, R) of radius R — [1/(^4(1+ | A l))''^^]- 

In the following, we also denote by f\ the complex analytic function defined for all 
A G on the disc D{0, R) by the Ao-series. 

Proposition 1.2.3 

For all 9, the series J2k>d^^^k{d), is an Ao-series with Ki — l/(d — 1), K2 — 



This is indeed a simple consequence of Bautin's method. For all 9, the norms of 
the polynomials A{9) and B{9) (seen as polynomials in (a, 6)), 



i/{d - 1), = i/W + = W + 1)] vd-i. 



Proof. 



A{9) = {xP + yQ)[cos{9), sin{9)], 



(1.14a) 



B{9) = {xQ - yP)[cos{9), sin{9)], 



(1.146) 



are bounded by: 



6 



I A{e) \< 2{d+l), I B{e) \< 2{d+l). (1.15) 

Write: 

dr/de = J2 r'^RkiO) = J](-l)^r'^+^('^-^)yl(^)S(^)^'. (1.16) 

k>d j>0 

Denote k = d + j{d — 1), then this yields: 

deg[Rk{e)] <l+j<l + [{k- d)/{d - 1)] < [{k - l)/{d - 1)]. (1.17) 

Furthermore, the norm of Rk{0) as polynomial in the parameters (a, h) is estimated 
as follows: 

I Rk{e) 1=1 A{e)B{ey \<\ A{e) \\ B{e) p< [2{d + i)f^-^^/^'^-^^\ (i.is) 

Theorem 1.2.4 

For all 6*, the series X]fc>d 3:'^t'fc(6'), is an Ag-series with K[ = l/{d—l), K'2 = 0, i^g = 
[2nKs/4Kl][K^-2C+{{K^-2Cf-Ki]^, 1/ K'^ = \K^-2C+{{K^-2Cf -Ki\/[2Kl], 
C = K^ + 2tt K^Kl. 

Proof. 

First observe that deg{vd{9)) = deg{Rd{d)) = 1 (The degree as polynomial in the 
parameters). Thus we have deg{vd{0)) < d/{d— 1). Assume inductively that: 

deg{v^{e)) < j/{d - l)forj = d, k - 1. (1.19) 
The recur rency relation (1.13) displays: 

deg[^;fc(^)] < max[deg{Bik[vd{(l>)),...,degivk-i{(l>)]Rii(l>))],i = 2,...,k. (1.20) 
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The equation (1-17) and the induction assumption yield: 



deg[B,k[vd{(l>),...,Vk-i{(f>)]Ri{(f>)]<K^{k-i) + K^i + K2<K^k. (1.21) 

This shows the first part of the theorem on the bound of the degrees of the coeffi- 
cients Vk{0)- For the second part of the proof related to the bound on the norms of the 
coefficients, we use standard methods and notations of majorant series. 

Definition 1.2.5 

The formal series ^{x) = X]a;>i ^k^^ with positive coefficients dominates the formal 
series ${x) — X^fc>i ^k^'' with positive coefficients if and only if for all coefficients 
^k <^k,k> 1. 

The series x + X]fc>2 | Vk{0) \ x'' is dominated by the series x + ^k>2 ^kx'' so that: 

k 

V^fc = 27r5^S,fc(V'd,->fc-i) \R^\■ (1.22) 

i=d 

Denote 

R{x) = J^KsKlx' = K3{K4xf/[l - K4X], (1.23) 

i>2 

then the proposition 1.2.3 yields: 

^ I i?i I a;*isdominatedby^(a;). (1-24) 

i>d 

The series ^{x) = x + X^fc>2 ^kx'^ is then dominated by the solution ^{x) solution 
of the equation: 

^(x) -x = 27r^[?^(a;)] = 2TTKz[K4^{x)f /[I - K#(a;)]. (1.25) 

At this point we have obtained that #(a;) is a solution of an algebraic equation of 
degree two. Estimates of the constant K'^ is then obtained by the distance to the first 
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zero of the discriminant and constant K'^ is then adjusted from the first term of the 
development. Easy computations give then the proof of the theorem. 

1.3 Definitions of the center set, of the Bautin ideal and of the Bautin 
index. Definition of ^-equivalence of series. 

Definition 1.3.1 

The center set C is the set of values of parameters (a, b) so that the corresponding 
vector field X has a center at the origin. 

Proposition 1.3.2 

The center set C is a real algebraic manifold. 
Proof 

The center set C is given as the zero set of the coefficients Vfc(27r) which are poly- 
nomials in the parameters (a, 6). 

From now on, it is appropriated to change of notations and denote L).(a, b) = Vfe(27r) 
the coefficients of the return mapping to emphasize their dependence in terms of the 
parameters (a, 6). We denote: 

r ^ L{r) = r + Ld(a, 6)/ + ... + Lfe(a, b)r^ + ... (1.26) 
the return mapping defined in (1.9) for 9 = 2tt. 

Definition 1.3.3 

The Bautin ideal is the ideal generated in the ring R[a, b] by the coefficients -^/.(a, b). 
Definition 1.3.4 
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The Bautin index is the first integer ko so that the polynomials Lrf(a, 6), LkQ{a, b) 
generate the Bautin ideal. 

The reader should be careful about the fact that the Bautin index does not depend 
only on the Bautin ideal. In its definition, we cannot substitute to the collection of 
the coefficients of the return mapping another system of generators of the Bautin ideal. 
Note that the existence of the Bautin index just follows from the fact that the ring 
R[a, b] is Noetherian. 

The local Hilbert's 16th problem is to find a bound depending only on d to the 
number of isolated periodic orbits of X in a neighborhood of the origin. Isolated periodic 
orbits of X defined in the domain of definition of the return mapping correspond exactly 
to the isolated solutions of the equation: 

L(r)-r = 0. (1.27) 

We gradually change into notations more pertinent to the general algebraic geom- 
etry setting. Let ${x, A) be an analytic series in x with polynomial coefficients in the 
parameter A = (Ai, A£)). 

Definition 1.3.5 

The Bautin ideal of the series <?(a;, A) is the ideal generated in the ring R[X] by the 
coefficients <?fc(a;, A). The center set of the series <P{x,X) is the zero set of its Bautin 
ideal. The Bautin index d of the series ^(x, A) is the minimal integer d such that the 
coefficients $i{x, A), ...$d{x, A) generate the Bautin ideal of the series ${x, A). 

Definition 1.3.6 

Two series ${x, A) and ^{x, A) with polynomial coefficients in the parameters A are 
said to be 0-equivalent if for all integers A; > 1, the polynomial <?fe(A) — ?f'fc(A) belongs 
to the ideal generated by (<?i(A), <?fc_i(A)). 
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It can be easily checked that this defines an equivalence relation. Also the definition 
yields the following: 

Lemma 1.3.7 

Two series ^(x, A) and W{x, A) with polynomial coefficients in the parameters A 
which are (/)-equivalent have the same Bautin index. 

1.4 The classical argument to show that the number of real zeros is bounded 

We recall here Bautin's argument to produce a bound to the number of real zeros. 
Let ${x. A) be an analytic series in x with polynomial coefficients in the parameters 
A and with Bautin index d. Assume that <5(0, A) = and ^{x, 0) = 0. 

Theorem 1.4.1 

There is a ball B e i?^, centered at 0, in the space and an interval / containing 
such that for all X E B, the number of zeros of <P{x, A) contained in / is less than or 
equal to d. 

Proof. 

Using the definition of the Bautin index d, we write: 

d 

^{x, X)^Y1 + ^ii^, (1.28) 

i=l 

with ^i{x, A) analytic in x, with polynomial coefficients in A such that ^i{0, A) = 
^i{x, 0) = 0. Assume that B and / are small enough so that (for instance) | ^i{x, A) |< 
1/2 on / X S. Then divide ^{x, A) by [1 + ^dix, A)] and rewrite: 



['P{x, A)]/[l + ^d{x, A)] = ^>i(A) + ^2(A)[1 + ^lix, X)]x^ + ... (1.29) 
Then from RoUe's lemma follows that the number of zeros of <P{x, A) is less than 
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1+ number of zeros of the derivative [^{x, X)]/[l + ^(i{x, A)]'. Write then this derivative 
as 

$2{X)[l + ^^^\x,X)]x + ...^d-i{X)[l+^P{x,X)]x'^-\ (1.30) 

Then repeat the process (sometimes referred to as the division-derivation algo- 
rithm). We obtain the result by an easy induction. 

It is clear that one cannot go far with this proof if we want an explicit control of 
the size of the domain (either of / or of B). For instance, we know that the coefficients 
<?j(A) are combinations of the d first coefficients: 

^i(A) = P(i,i)(A)^i(A) + ... + P(,,d)(A)^d(A), (1.31) 

but how to control the size of the "quotients" P(i fc)(A) which enter in the construction 
of the ^i{x, A)? We will see how to bypass these difficulties in the next chapter using 
complex analytic methods such as the polynomial Hironaka division theorem. 

1.5 Formulation of the problem in terms of projection of analytic sets. 

Let <P : Rx i-^ i? be an analytic series with polynomial coefficients: 

^{x, \)^x + $2{X)x^ + ... + ^fe(A)a;'= + ... (1.32) 
We consider the subset E C R x R^ defined as the zero-set of <P{x, X) — x: 

E=[{x,X)/^x,X)-x = 0]. (1.33) 

Let TT : U R^ be the restriction to U of the natural projection tt : Rx R^ R^. 
The center set associated to the analytic series ${x, X) — x defined in 1.3.5 is the set 
C C R^ of parameters A such that the fiber of the projection 7r~^(A) is contained in 
the set E. The theorem L4.1 displays the following geometric interpretation: 
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Theorem 1.5.1 



There is a neighborhood I x B oi (0, 0) in x such that for all points A of 5 
the number of isolated points of the fibers 7r~-'^(A) restricted to E is less than the Bautin 
index d of the analytic series ${x, X) — x. 

1.6 A S'^-action on the space of parameters and a system of S"^- invariant 
polynomials which generate the Bautin ideal. 

We can change the coordinates in the plane {x, y) ^ {x\ y') under the action of the 
rotation group. This leaves invariant the linear part of the vector field X and thus it 
induces an action of the rotation group on the space of parameters (a, h). This linear 
action induces an action on the ring of polynomials R[a, h] on the space of parameters. 
A polynomial is said to be invariant if it is fixed under this (S'^-action. The center set is 
obviously invariant under the S'-'^-action because the fact, that the vector field X has all 
its orbits periodic in a neighborhood of the origin, does not depend of the choice of the 
system of coordinates on this neighborhood. It is natural to ask if the coefficients of the 
return mapping are /S^-invariant polynomials. The answer is no, but it is possible to 
find another system of generators of the Bautin ideal which are (S'^-invariant and so that 
the associated analytic series is 0-equivalent to the return mapping. This was proved by 
([6]) using the approach of the successive derivatives. For d = 2, 3 it was proved by H. 
Zoladek ([15])using the classical Bautin approach. We recall here his proof and derive 
the result in its full generality for any d, (d > 2). 

Let us consider the ring of polynomials in (a, 6), [^, (sin(^), cos(^)]: 

i?[a, h] (g) R[e, (sin(^) , cos(^)] . 

The rotation group acts on the ring i?[a, 6] as displayed previously and so it acts on 
i?[a,6] (sin(^), cos(^)] in an obvious way (just extend the action by the identity 

on the second factor). We denote p i— > * p this action where belongs to and p 
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belongs to R[a,b] ^R[d, (sin(^), cos(^)]. Now 5'-^ acts also on the ring 

R[9, (sin(^),cos(^)] 

in a natural way: 

* (6») = 6* + 0, * (sin(6'), 008(6*)) = (sin(6' + 0), 008(6* + 0)). (1.34) 
This action can be extended to an action p 1— > (f)*p on the ring 

R[a, b] (^R[e, {sm{e), cos(6')] 

by the identity on the first factor. 

Definition 1.6.1 

An element p of the ring R[a, b] I^R[9, (8in(^), cos{9)] is said to be covariant if for 
all 0: (f)* p — (f)*p. 

Lemma 1.6.2 

Product and sum of two covariant polynomials are covariant polynomials. 
The proof of this lemma is very easy and is omitted. 
Lemma 1.6.3 

The coefficients Rk(0) of the equation (1.8) are covariant polynomials. 
Proof. 

The trigonometric polynomials A{9) and B{6) are clearly covariant polynomials 
as seen from the definition of the actions. Now from Lemma L6.2 it follows that the 
products A{d)B{6)^ are covariant polynomials. 
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It is now necessary to be more careful with the inductive construction of the coef- 
ficients Vk{0). We initiate the discussion with the term of lowest degree Vd{0) defined 
by the condition: 

vM = Rd{0)- (1.35) 
The term RaiO) to be integrated is a trigonometric polynomial. Thus this yields: 

Vd{e) = zd{e) + Sd{9), (1.36) 

where Zd is a constant and Sd{P) is a trigonometric polynomial with Sd(0) = •Sd(27r) = 0. 
Note that 

/•27r 

Zd = (l/27r) / Rd{(l>)dcl>, (1.37) 
and the first coefficient of the return mapping is: 

Vd{27:)^2nzd^ / Rd{<j))d(p. (1.38) 
Jo 

Next step of the recurrence: 

v'd+iiO) ^J2^^d+i[vd{0)mO), (1-39) 

i=d 

leads to: 

Vd+i{0) = Zd+iO + ri2i^2 ^ sd+i{e), (1.40) 
where Sd_|_i(6') is a trigonometric polynomial so that: Sd_)_i(0) = Sd+i(27r) = 0, 

^d+i = (l/27r) / [V5id+i[sd(0)]i?i(0)]d0, (1.41) 
-^0 i^d 

and where the coefficient Rd+i is proportional to the coefficient Vd+i- 
The general step k of the recurrence displays: 
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where Sk(0) is a trigonometric polynomial so that: Sfc(O) = Sfe(27r) 
such that: 



(1.42) 

= 0, 



k 



(1.43) 



.27r k 



Zk = (l/27r) / [^5,fc[sd(</'),...,Sfc-i(</>)]i?z(</>)]#, 



(1.44) 



and where the coefficients , j = 2, — (i+ 1 belongs to the ideal generated by 
the preceding coefficients za-, ...,Zk-i. The construction of the coefficients Zk yields: 

Lemma 1.6.4 

The two series L{x) — x — '^k>aVk{'2TT)x^ and ${x) — x — Ylik>d'^kX^ are 0- 
equivalent. 

We prove now that the polynomials Vk are ^^-invariant. The proof splits into 
several lemmas of independent interest. 

Lemma 1.6.5 

Let Tk{9) be an element of the ring R[a, h] R[d, (sin(^), cos(^)] which is a covariant 
trigonometric polynomial, then the polynomial defined by: 




(1.45) 



is an invariant polynomial. 



Proof. 



Let 6 be an element of ^ consider 6 * (v^): 
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/■27r 

e*{vk)^ / e*Tk{(t))d(t>. (1.46) 

Jo 



The covariance property of yields: 



/»27r nliz /•27r /•27r 

/ e*Tk{(i>)d(j)^ e*^Tk{(j))d(i) ^ n{ct> + e)d(t>^ Tk{(t>)d(t). (1.47) 

7o JO JO Jo 

Lemma 1.6.6 

Let Tk{9) be a trigonometric polynomial without constant term which is covariant. 
The only trigonometric polynomial Sk{0) without constant term so that: 



si{e) = n{e\ (i.48) 

is covariant. 
Proof. 

Write 

Tk{d) = Y,Tki{a,h)e'''. (1.49) 

Let us consider the first action of S^. Let (j) be an element of S^. Under the 

action of this element, each term Tki(a,b) gets multiplied by the factor e**^. Now the 
polynomial Sk{0) writes: 

Skie) = J2[Tkiia,b)/iy''. (1.50) 
I 

Its covariance follows easily from this expression. 



We conclude this paragraph with the theorem 
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Theorem 1.6.7 



The displacement function L{x) — x = ^fc>(i t'fe(27r)x'^ is 0-equivalent to a series 
>d^kX with S'^-invariant coefficients Zk- 

Proof. 

The polynomial RdiO) is covariant (cf. Lemma 1.6.3) and thus the coefficient Zd is 
invariant as a consequence of formula (1.44) and of the Lemma L6.5. The polynomial 
Sd is covariant as a result of the Lemma L6.6. Assume inductively that the polynomials 
SdO, ...,Sk-\{9) are covariant. Then Lemma L6.3 and Lemma L6.6 and the formula 
(1.43) yield the fact that Sk is covariant and the Lemma 1.6.5 imply that the coefficient 
Zk is invariant. 

II-A GLOBAL AND COMPLEX GENERALIZATION OF BAUTIN'S 
THEOREM. 

Bautin's method as displayed in the preceding chapter does not allow to produce a 
precise estimate on the size of the domain on which the number of zeros of an Ao-series 
(such as the displacement function) is controlled. It is not surprising that to get such 
an estimate, complex analytic methods are more appropriated. We begin this chapter 
with general considerations on Bernstein classes. 

II. 1 Finiteness properties, Bernstein inequality. 

In this part, the definitions of the Bernstein classes Bi and B2 are displayed and 
their equivalence is proved. Then it is shown that the number of zeros of functions of 
these classes is explicitly bounded. The closed disc centered at e C of radius R is 
denoted Dn. This paragraph relies entirely on ([8]) and ([9]). 

Definition II. 1.1 
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Let i?>0,0<Q!<l and K > and / holomorphic in a neighborhood of Dji. 
The function / belongs to the Bernstein class B\^^ j^ if and only if: 

max[| f{z) I, z e Dii]/max[| f{z) \,ze Do^r] < K. (2.1) 

Definition II.1.2 

Let N be an integer, i? > and c > and f{z) = X]j>o O'i^^ be an analytic function 
on a neighborhood of e C. The function / belongs to the Bernstein class Bf^ r c 
and only if for all j > N: 



aj 



W < cmaxi=o,...,iv(| | R'). (2.2) 



Theorem II. 1.3 



Let / be an element of ^ ^. Then / is analytic on the open disc Dr, and for 
aU R' < R ond a < 1, and K = a-^[l + a{l - a^)/{l - a) + c/3/(l - /?)] where 
P = R'/R, f belongs to Bj^, ^ ^. Conversely if / belongs to 5)j ^ ^, then / e Bf^ ^ ^ 
with N = [(log2(K) - log2(l -a) + l)/(log2(l/«))], c = [K{2K + 1)/(1 - a)^. 

Proof 

If / e Bjf ^ the convergence of the series f{z) = J2i>o ^i^^ ^^e disc Dr is 
consequence of the Cauchy-Hadamard formula and of the inequality: 

I an \< cmaxi=o,...,iv(| a, | R'){R-''), (2.3) 

for large n. Denote m = max(| f{z) \,z & D^r, Cauchy's inequalities yield: 

\ai\<m/{aR')\ (2.4) 

For i = 0, N, this displays: 

\ai\R' < m/{aR'/Ry < [m/ {R' / R)^] (2.5) 
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and then for all j > N + 1: 



I aj I < c[m/a^{R'/R)^]. (2.6) 

This yields: 

N 

i\ f{z) \,zeDR,<J2 + E l^^l^"' (2.7) 

i=0 j>N+l 



max 



TV 

<mE(l/ai?')*i?'^ + c[m/a^(i?7i?)^] ^ {R'/Ry (2.8) 

i=0 j>N+l 

, and thus: 

max(| f{z) \,ze Dr.) < [m/(«^)][l + «(1 - «^)/(l - a) + cp/{l - /?)] (2.9) 

as expected. Conversely, assume that / belongs to Bj^^ j^, f{z) = X]i>o'^«'^*- 

N be an integer and a = maxi=o,...,Ar(|ai|i?*) and b = max(|i?(z)|, z G D^r), where 

-R(^) = Y^i>N+i ^i^^- Consider Pn{z) = Ya=q ^i^' and: 

N 

max(|Pjv(-2)|,2;eDaii) <E|a,|(ai?)^ (2.10) 

1=0 

max(|PAr(2;)|,2; e D^r) < a /{I - a). (2.11) 

This yields: 

uia;^{\f{z)lz e D^r) < K[h+ {a/{l + a)]. (2.12) 
Cauchy inequalities yield: 

\aj\R' < K[h + [a / (1 + a)]. (2.13) 
The following bound for h may be displayed: 

h< J2 [Wj\{(^Ry <K[b+{a/{l + a)][a''+^/{l-a)], (2.14) 

j>N+l 

b[l - (Ka^+V(l - «)] < Ka[a^+^/{l - af]. (2.15) 
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Fix N so that: 

(Ka^+V(l - «) < 1/2; (2.16) 

this is achieved for instance if: N = [(log2(ii') — log2(l — a) + l)/(log2(l/Q;))]. Such a 
choice for N yields: 

b < 2Ka[a^+^/{l - af] < 2Ka/{l - af. (2.17) 

Now this yields: 

\aj\W < K[2aK/{\ - af] + [(7/(1 - a)], (2.18) 

\aj\R> <K(2K + l)/{l-af; (2.19) 
the theorem now follows with: c = [K{2K + 1)/(1 — ct)^. 

The number of zeros of a function which belongs to a Bernstein class is bounded 
accordingly to what follows. 

Theorem II. 1.4 

If / G Bj^ ^ j^, then the number of zeros of / in the disc D^n is less than: 

log2(K)/[log2[(l + a2)/2a]] (2.20) 

Proof: 

Assume that / displays n zeros ^i, ^„ in the disc D^iji. Denote g the holomorphic 
function on the disc Dji defined as: 

n 

9{z) = f{z) H [(i?2 - zzk)/R{z - Zk)]. (2.21) 

k=l 

The maximum principle yields: 

n 

max(|^(^)|,^ e Dr) = max(| f{Re'') J][(i?2 _ Re^^ Zk) / R{Re'' - Zk)] \,e e [0,27r]), 

k=l 

(2.22) 

21 



= max(|/(i?e^^) || 1[[{R^ - Re'' Zk) / R{R - e'' Zk)] |, ^ G [0, 27r]), (2.23) 
fc=i 

= inax{\f{z)\,zeDR). (2.24) 

Consider now: 

n 

max(|^(^)|, z e Dan) = max(| f{aRe"^)) \[[{R^ - aRe'^ Zk) / RiaRe'^ - Zk)] |) (2.25) 

k=l 

n 

> max{\f{z)\,z e D^r) JJ min[| R^ - aRe'^ Zk) / R{aRe'^ - Zk) \ ■ (2.26) 

k=l 

Write now separately each quantities: 

I (i?2 - aRe'^Zk)/RiaRe'^ - Zk)] \ (2.27) 
and Zk — RkG^^'' ■ The minimum of the quantity when 9 — 9k varies in [0, 2tt] is 

{R + aRk)/{aR + Rk). (2.28) 

The minimum of this quantity when < Rk < cxR is 

{l + a^)/2a. (2.29) 

This yields the following inequality: 

max{\g{z)lze DaR) > max(|/(^)|, ^ e Dafl)[(l + «2)/2«]^. (2.30) 
The inequalities: 

max(|/(^)|,z e Dr) = majc{\g{z)\,ze Dr) > max(|^(^)|, ^ e Dau), (2.31) 
display: 

K >[{l + a^)/2a]'', (2.32) 
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and thus: 



n<log2(K)/[log2[(l + a2)/2a]]. 



(2.33) 



Lemma II. 1.5 

Let / e S^^c and R" = i?/[23^max(c, 2)], the number of zeros of / in the disc 
Dn" is less than N. 

Proof. 

Set R' = i?/[max(c,2)] and a = 2"^^. This displays: 

[(1 + a')/2a] = (2^^ + l)/{2^^+^) > 2^^-\ (2.34) 

K < 23^-^+2. 
This yields that the number of zeros of / is less than: 

(3N^ + 2)/(3N -1) < N +1, (2.35) 

as soon as N > 2. 

II. 2 The Hironaka polynomial division theorem 

The precise statement of the following theorem was provided by P. Milman. The 
terminology "Hironaka polynomial division theorem" is also due to him. This technical 
result allows to clarify the presentation of the "Quantitative Bautin's theorem" first 
proved in ([8]). 

Theorem II. 2.1 

Let / be an ideal of C[Ai, Ad]. There is a system of generators gi, Qs of the 
ideal / and constants C and Ci such that for all elements / of degree k of /, there is a 
decomposition: 
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with 



s 



(2.36) 



deg{4>^) < deg{f) = k, (2.37a) 

and 

|<^i|<CCi"|/|. (2.376) 

We include here a fuU proof of this theorem. The non classical part of this division 
theorem is displayed in the control of the norms (equation 2.37b) of the quotients. The 
classical proof uses the inversion of an operator defined on Banach spaces of analytic 
functions but we do not follow these lines here. 

Definition II.2.2 

In the following, a total ordering < on is said to be compatible with the addition 

if: 

i) For all indices a e N^, (3 e , then a < a + {3, ii) For all indices a^, Q;^,/3, 
+ j3 < -\- j3 ii and only if < . 

Here, we choose (for instance) the total ordering defined in such manner: 

Q! < /? if and only if a.\ + ... + ckd =| a |< /?i + ... + =| /? | or if | a | = | /5 |, 

then C{oi) < C(/3) where C{a) = X^i^o*^*^*' ^ is a transcendent number, < e < 1. 

Note firstly that < defines a total ordering compatible with the addition. Indeed, the 

transcendent nature of e yields 

C{a) = C{(3)^a = /3 (2.38) 

Definition II.2.3 
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Let / e C[A], / ^ 0, f = fa^°'- The largest exponent a so that fa^Ois called 
the privileged exponent of / and is denoted exp{f). The monomial [/q,A'*,Q! = exp{f)] 
is called the initial monomial and denoted In{f). 

Definition II.2.4 

Given s polynomials gi, Qs, the associated partition of A^^ is defined as follows: 
A^ = exp{gi) + N^, A, = expig^) + - [j.^, A^, A = - Ul=i A. 

Definition II. 2. 5 

Let / be an ideal of C[Ai, A^]- Consider the set: 

exp{l)^{exp{f)jel). (2.39) 

It can be shown that this set has finitely many extremal points: 

a^,...,a'. (2.40) 

Choose gi, gs in the ideal / so that exp{gi) = a*, z = 1, ...s. Such a set of polynomials 
is called a standard basis, Hironaka basis or Grobner basis (of the ideal / relatively to 
the ordering <). 

Proposition II. 2. 6 

Let / G C[A], deg{f) = k, and gi, gs € C[X\, there is a constant C which depends 
only on the polynomials gi and there are unique hi, /ig, h e C[A] such that: 



/ = higi + ... + hsga + h, (2.41a) 

hi = ^ hi,a>^'^ ^ a + exp{gi) e A, (2.416) 

h^Y^ /laA" ^ a e (2.41c) 

Max{\ hi\,\h\)< CC^ I / I, (2.41d) 
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Max[deg{hi), deg{h)] < k. 



(2.41e) 



Proof: 

Given /, denote In{f) = /aoA"o. If ao e ^ set /i^^) = /aoA"" and /if ^ = 0, i = 
1, s, /(I) = / - This yields | /(i) |<| / | , deg{f^^^) < deg{f) and | /i^^) |<| / | 

may belong to several Aj, choose one). Then set: 
hf^ = Oi{jj^i and h^^^ = 0. Write ao = exp{gi)+Po, In{gi) = 7iA«^P(»*), then denote 
hl'^ = /«oA^V7i and /(i) = / - h^^g^. This yields: 

\h[^^\<C\f\ (2.42) 

where: 

C = max(l/7,),j = l,...,s (2.43) 

deg(/if ^) ^\Po\<k- min(| exj9(^j) |) (2.44) 

and thus this displays: 

|/^'M<I/I+CG'|/|<(1 + CG)|/| (2.45) 

where: 

G = max(| gi \),i = 1, ...,s. 

Note furthermore that: 

deg(/(i)) < max[deg(/), deg(/if ^^i)] (2.46) 

and 

deg{h^^g^) =1 Po I +deg(^i) <| /3o I + I ex^^^) |<| ao |< k (2.47) 

and thus 

deg(/(^)) < k (2.48). 

Then, repeat the whole process with f^^^ in place of /. Note that the privileged exponent 
of /(^^ is strictly less than the privileged exponent of /. The process stops ultimately 
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when the privileged exponent becomes zero. The choice of the ordering yields that 
the number of steps involved is less than k{l + e^~^). This yields the result of the 
proposition with: 

Ci = {1 + 00)^+'''" (2.49) 

The theorem follows of the preceding proposition and of the property that if a 
polynomial / belongs to an ideal / and if gfi, fif^ is a Grobner basis of the ideal /, then 
the division of / yields h = 0. (See for instance Lejeune-Jalabert ([10])). 

II. 3 The main theorem 

Theorem II.3.1 

Let fx{x) = 5^fe>i /fc(A)a;'^ be an ^o-series. For any A e C^, the function fx 
belongs to the Bernstein class ^ ^ where d is the Bautin index of the series f\{x), 
R = [(ClA)^l^^4]"^ c = [MCK3(CiX)^^]/R'^, if i? < 1, c = [MCK3(CiX)^^]/R, if 
> 1, A = max{l,\ A |). The constants Ki,...,K4 are those which appear in the 
definition (1.2.1) of an Ag-series. The constants C and Ci appears in the "Hironaka 
polynomial division theorem" and M is the norm of the matrix of the change of basis 
(fk) to a Grobner basis (gi) of the Bautin ideal of the series fx{x). 

Proof 

Write first the condition for f\{x) to be an ^o-series as follows: 

deg[fk{X)]<Kik + K2, 
\fk\<K^Kl 

Denote / the ideal generated by the d first coefficients /i(A), /d(A) and write: 

s 

fkW = J2^4>';W9jW]. (2.50) 
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The Hironaka polynomial division theorem yields: 

I 0,^(A) |< I I x^K^k+K,) (2.51a) 

< [CK3C^'][C^'K^]''X^^'''+^^\ (2.516) 

The relation: 

d 

9j{X) = J2^ji{X)MX), (2.52) 
1=1 

yields: 

I MX) \< [CKsCf ^A^^][C^i^4A^^]'= J2 I 9jW I, (2.52) 
< M[CK3C^A^^][C^i^4A^^]'=niax[| /i(A) |, | /d(A) |]. (2.53) 

Write now: 

i?=[(CiA)^iK4]"\ (2.54) 

the equation (2.53) yields: 

I fk{X) I R'^ < [MCK3(QA)^V^1niax[| /i(A) |, | fa{X) \]R''. (2.55) 
li R < 1, this displays: 

I fk{X) I i?'^ < [MCi^3(C^iA)^V^Vax[| /i(A) I R, | /^(A) | i?'^]. (2.56) 
If i? > 1, the equation (2.53) yields: 

I /fc(A) I i?'^ < [MCK3(CiA)^V^]niax[| /i(A) |, | fd{X) \]R, (2.57) 

and thus: 

I fk{X) I R'' < [MCK3(CiA)^V^]niax[| /i(A) | R, | fd{X) \ R% (2.58) 

This concludes the proof. 

The theorem II. 3.1 and the lemma II. 1.5 now implies the following: 
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Theorem II.3.2 



The number of zeros of fx{x) in the disc Dji" is less than d — 1 with 

R" = i?/[23('^-i)max(c, 2)]. (2.59) 

III. Extension to any dimension of Bautin's theory 

These last years, the dynamics of plane systems was extensively studied and sev- 
eral new techniques were developed. Some are specific to 2-dimensional systems but 
often these methods can be appropriately extended to multi-dimensional systems. The 
algorithm of the successive derivatives was derived some years ago ( [5] ) to find the first 
non-vanishing derivative (relatively to the parameter e) of the return mapping (near the 
origin) of a plane vector field Xq + eXi of type: 

d 

Xo + eXi = xd/dy-yd/dx + e ^ [aijxyO/dx + bijxyd/dy]. (3.1) 

The algorithm was then used in the center- focus problem (cf. [6]), which directly 
relates to Hopf bifurcations of higher order and to several other problems on limit cycles 
of plane vector fields. How to extend appropriately this situation in any dimension? 
We have to perturb a dynamics which is integrable and displays only periodic orbits. 
Assuming that the perturbation depends of finitely many parameters (say for instance it 
is polynomial), we expect also that the perturbed system displays a first return- mapping 
which is analytic with a Taylor expansion with coeflBcients which depend polynomially 
of the parameters. This return- mapping should label all the periodic orbits (at first 
return) of the perturbed system by its fixed points. The principal aim of this paragraph 
is to present a framework where such demands are realized. In this framework, a 
generalization of the algorithm of the successive derivatives is provided. 

III-l Controlled Nambu dynamics and (*)-property. 
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Let / = (/i, ...,/n-i) : — > R'^~^ be a generic submersion (meaning that / is 
a submersion outside a critical set f~^{C), where C is a set of isolated points). Let 
Q — dxiAdx2A...dxn be a volume form on i?". Consider the vector field Xq such that: 

iXodxiAdx2/\---/\dxn = dfiA...Adfn-i- (3-2) 
The functions /j, (i = 1, n — 1) are first integrals of the vector field Xq: 

dfif\iXodxif\dx2f\---f\dxn — {Xo.fi)dxiAdx2/\---/\dxn — dfiAdfiAdf2A...Adfn-i — 0. 

(3.3) 

This type of dynamics is well-known in Physics and named Nambu's dynamics 
([12]). 

For c varying in a neighborhood of 0, assume that the curves f~^{c) have a compact 
connected component 7c. Let U he a, small neighborhood of the zero-section of the 
normal bundle to 70. For c small enough, the curves 7c are closed periodic orbits of Xq 
and they cut transversely E. Choose c as a coordinate on the transverse section E to 
the flow of Xq. Lastly, assume that there are 2-forms oji such that: 

tXo'^i = 4fi;^ = 1, 1- (3.4) 

There are of course different possible choice of the forms oji and accordingly different 
possible perturbations have to be considered. If the condition (3.4) is fulfilled, we will 
say that the singularity of the Nambu dynamics (3.3) is controlled (or alternatively that 
the Nambu dynamics itself is controlled). 

The appropriated extension of the (*)-property first discussed in ([5]) is presented 
in the following. 

Definition III. 1.1 

Let / = (/i, ...,fn-i) : —>■ R^~^ be a generic submersion. Assume that /~^(c) 
contains a compact curve 7c. The application displays the (*)-property if for all poly- 
nomial 1-forms uj such that 
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a; = 0, 



(3.5) 



for all c; there exist polynomial such that: 



<^ = 9\dii + ... + fifn-l 



dfn-1 + Ci-R. 



(3.6) 



It was proved in ([5]) that the function /i : i?^ — > R^ifi : {xi,X2) — > {xf + x^) 
displays the (*)-property. Several generalizations were proposed after but the core of the 
argument in the computation of the successive derivatives is captured in this notion. The 
generalization proposed in this article provides a new presentation of the (*)-property 
which seems interesting as well for the 2-dimensional case. Indeed, the definition of the 
vector field Xq given in the preceding introduction yields the: 

Proposition III. 1.2 

Let a; be a 1-form such that u!{Xo) = 0, then there are functions gi, gn-i so 
that: 



Note that the condition u{Xo) = 0, equivalent to u)/\dfiA...Adfn-i = 0, yields 
LV = gidfi + ... + gn-idfn-i where the coefficients g^ are obtained as ratio of minors of 
the Jacobian matrix of the fj. 

This displays an alternative to the (*)-property now presented as follows: 

Proposition III. 1.3 

A generic submersion / : i?"^ — > R'^~^ displays the (*)-property if for any polyno- 
mial 1-form uj such that 



= gidfi + ... + gn-i 



dfn-i- 



(3.7) 




(3.5) 



31 



for all c; then there exists a polynomial R such that: 



u;{Xq) = Xo.R. (3.8) 

Such a function R can be (in principle) constructed with the following pattern. 
Choose R arbitrarily on the transverse section E, then extend R to the whole tubular 
neighborhood of 70 saturated by the orbits 7c by integration of the 1-form u along the 
orbits of Xq. 

III-2 The successive derivatives of the first return mapping of the per- 
turbed system. 

Now perturb Xq into X^ = Xg + eXi. Let M be a point of S close to and let 7e be 
the trajectory of X^ passing by the point M. The next first intersection point of 7e with 
E defines the so-called first return mapping of X^ relatively to the transverse section 
U: c I — > L{c,e). The mapping L is analytic and it displays a Taylor development (in 

L{c, e) = c + eLi(c) + ... + e'=Lfc(c) + 0{e)''+\ (3.9) 

The expression of the first coefficient Li{c) is classical and belongs to the lore of 
bifurcation theory. With the vector field and the 1-forms u)i (cf. [F]), introduce the 
1-forms: 

tXe^i = i^Xoi^i + ^i^Xii^i = dfi + eux^tOi- (3.10) 

Definition III.2.1 

The perturbation X^ of the controlled Nambu dynamics is said to be admissible if 
for all the 2-forms Wj, the 1-forms ix-^^i have polynomial coefficients. Note that has 
said above this displays different admissible perturbations depending of the choice of 
the forms a;^. 

32 



Recall that the parameter c chosen as coordinates on the transverse section S is 
the restriction of the functions / = (/i, fn-i) to the section. 
Then the i^'^-component of Li{c) is equal to: 

Li,i{c) = / ixi<^i- (3-11) 

-^70 

Assume now that the first derivative Li{c) vanishes identically and that the sub- 
mersion / displays the (*)-property then there exist gij and Ri such that: 

iXi<^i = ^9ijdfj + dRi. (3.12) 

3 

Following the lines of the algorithm of the successive derivatives, the expression 
(3.12) yields: 

i^2,i(c) = - / ^gij^x.ujj. (3.13) 

This is indeed the second step of a general recursive scheme which displays as 
follows: 

Assume that all the /c^'^-first derivatives of the return mapping of the perturbed 
vector field vanish identically. This yields: 

LkAc)- [ 5^4-^x,a;, =0. (3.14) 
The (*)-property yields new functions gfj,R^ such that: 

Y.g'^rhx.ujj = Y.g%df^ + dRK (3.15) 

3 3 

This yields the following expression of the {k-\- l)*''-derivative of the return mapping of 
the perturbation: 

Lk+i,i{c)^ ^g'^jix.ujj. (3.16) 
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The algorithm imphes of course the first 



Theorem III.2.2 

Let Xq be a controlled Nambu dynamics which displays the (*)-property and let Xi 
be an admissible perturbation. Then the perturbed dynamics has an analytic first 
return map. The coefficients of the Taylor expansion of this return mapping depend 
polynomially on the coefficients of the perturbation. 

From the general theory of projections of analytic sets ([11]), it now follows: 

Theorem III.2.3 

There exists a uniform bound to the number of isolated periodic orbits, which 
correspond to fixed points of the first return mapping of Xq + eXi which intersect the 
transverse section E in the neighborhood of 0. 

The general framework presented here should of course be illustrated with specific 
examples (of dimension larger than 2). Some have been worked out recently by Seok 
Hur (Paris VI) and will be matter to further publications. 
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